Transverse localization and slow propagation of light 
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The effect of finite control beam on the transverse spatial profile of the slow light propagation in 
an electromagnetically induced transparency medium is studied. We arrive at a general criterion in 
terms of eigenequation, and demonstrate the existence of a set of localized, stationary transverse 
modes for the negative detuning of the probe signal field. Each of these diffraction-free transverse 
modes has its own characteristic group velocity, smaller than the conventional theoretical result 
without considering the transverse spatial effect. 
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Ultraslow propagation of light fields has been an active 
research field recently. Controlling the light propagation 
in atomic and solid-state media is important in both the 
fundamental theory and practical applications of nonlin- 
ear optics 1\. Slow light propagation experiments have 
been reported by using ultracold atoms 0,0, 0, hot 
atoms 0, El, rare-earth ion doped crystal pj, ruby 0], 
and alexandrite crystals 0]. The use of electromagnet- 
ically induced transparency (EIT) to obtain slow light 
propagation is one of the most important techniques. In 
an EIT medium when a control laser is applied to 
an appropriate transition, a weak probe signal pulse may 
have small absorption and steep dispersion. Due to this 
steep dispersion, the group velocity of the signal pulse 
can be reduced to several orders smaller than the light 
speed in vacuum 0, 0, H Further, by chang- 
ing the intensity of the control laser, it is possible to 
reversibly stop the signal pulse 0. Stop light propa- 
gation has been observed in cold and hot alkali vapors 
|12L Il3| . Possible applications of slow light include the 
enhancement of nonlinearity entanglement of 

atomic ensembles or photons [171 1 1 8| , quantum memories 
[l9|. and optical information processing pfj|. 

Almost all theoretical treatments so far involve the as- 
sumption of an effective infinite transverse spatial varia- 
tion of control field, and how a finite transverse profile af- 
fects the light propagation in an EIT medium is yet to be 
addressed. An exception is the work in Ref. 0. The au- 
thors investigated the transverse localization of the sta- 
tionary probe pulses with a pair of counter-propagating 
control fields, and showed it is possible to realize a three- 
dimensional confinement of the probe pulses |2l|. The 
focus of Ref. |2l| was how to localize light pulses in three 
dimensions. In this paper, we are interested in the trans- 
verse effects on slow light propagation. In particular, 
we arrive at a general criterion in terms of eigenequa- 
tion for the transverse stationary modes of the signal 
field. These transverse modes exist only at negative de- 
tuning frequency region, and do not undergo diffraction, 
thus keeping their profiles during the propagation. The 
group velocity of each stationary mode will also be stud- 



ied quantitatively. 
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FIG. 1: The three-level A-type EIT medium. Q. is the Rabi 
frequency of the classical control field, and E is the electric 
field of the probe signal pulse. 

We start with the three-level A-type atoms shown in 
Fig-ffl The medium of length L consists of an ensemble of 
N atoms. The ground state |0) and the metastable Stokes 
state |1) are coupled individually with the excited state 
1 2} via a weak probe signal pulse and an intense control 
laser field, respectively. The latter is treated classically 
and assumes the form £ c (z,t,r) = Q(z,t,r)e l ( kcZ ~^ ct \ 
where Q(z,t, r) denotes the Rabi frequency and k c = 
lu c /c, with the carrier frequency, to c — u>i2, being set to 
be on resonance with the Stokes transition. The weak 
signal field shall be treated as a quantum field : 



£(z,t,r) = E(z,t,r)e i(koZ - iW02t \ 



(1) 



where ko = lo^/c. The probe carrier frequency oj s is as- 
sumed to be close to the |0) — ► |2) transition frequency 
cjo2- In Eq. QJ, E(z,t,r) denotes the slowly varying sig- 
nal field envelop operator. Following Ref. till but includ- 
ing now also the finite control beam effects, the signal 
pulse in the paraxial approximation and slowing varying 
amplitude approximation can be describe by the propa- 
gation equation, 
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E = igNa 2(z,t,r). 
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Here. 



V\ = V 2 - 



d 2 



is the transverse Laplacian, <? = 
/i02 J 2heaV ' with M02 being the transition dipole moment 

and V the quantization volume, denotes the atom-held 
coupling constant for |0) — > |2), while 002 is a slowly 
varying collective operator of the atoms. In general, 

&ij(z,t,r) — jpJ2^=i r N)Kjl e ~ JWii 'j where the sum runs 
over the effective number of atoms in a small but macro- 
scopic volume around position r |2l| . 

We shall be interested in the case in which the Rabi 
frequency of the signal held is much smaller than that 
of the control field Q and the number of input probe 
photons is much less than that of atoms. In the adiabatic 

approximation, &02(z, t, r) = ^§|j— EHHH l; leading 
Eq. © to 
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In most experiments, the control field is continuous and 
has a cylindrical symmetry transverse spatial profile f2(r) 
that changes little in the propagation direction. To study 
the transverse effects of the probe signal field in this case, 
we consider the expectation value of E(z,t,r) in terms 
of 



E(z, t,r) = ip(r)e 



im9 ^i(Az- 
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(4) 



with the quantum number m = 0, ±1, ■ • • for the orbital 
angular momentum of the signal field The signal 

wavevector mismatch A = k s — fco along the z-direction 
will be determined as the function of frequency detuning 
S = uj s — CJ02 • Equation © can then be reduced to 



dr 2 



1 d 2k g 2 NS 
r dr cfl 2 (r) 



^W=WW, (5) 



where the eigenvalue (3 = 2ko(8/c— A) defines the disper- 
sion relation between A and 8 for each transverse mode 
of the signal field. The physical boundary conditions for 
Eq. © are ^1(00) = V'(0)| m=0 = V>(0)|„^ - 0. 

It is noted that Eq. © may be considered as a prop- 
agation version of Eq.(2) in Ref. |^] ( by s etting f2_ = 
there). However, our work and Ref. |2l| have different 
physical background. Ref. was focused on how to re- 
alize light localization. The novel result of this paper is 
Eq. JSJ. It unambiguously shows the general properties 
of a complete set of transverse invariant eigenmodes for 
slow light propagation in the EIT medium. 

For a negative detuning (8 < 0), the eigenequation 
(|5|) for each integer value of m is quantized. The re- 
sulting transverse modes and eigenvalues are denoted as 
ipmn(r) and f3 mn , respectively; with n — 1,2, ••• and 
< Pmi < Pmi < • ■ ■ for each m = 0, ±1, • • • , ±(n — 1). 
Physically, the negative frequency detuning leads to the 
decrease of the refractive index from the optical axis and 
the production of an effective waveguide for the signal 
field. The group velocity for each stationary transverse 



mode of the signal field can be evaluated via its eigen- 
value as 
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For a given control field with an arbitrary transverse pro- 
file and the relevant parameters for the optical medium, 
one can solve Eq. (JSJ) to obtain all eigenvalues (3 mn and 
the corresponding transverse modes Tpmnir). Any local- 
ized input probe signal pulse E(0,t,r) can be expanded 
as the superposition of these transverse modes; each of 
these modes then propagates at its own distinct velocity, 
and thus, the signal pulse of superposition changes its 
spatial-temporal profile during the propagation. Clearly, 
if the input pulse is characterized by a single transverse 
mode it will not undergo diffraction and remains the ini- 
tial shape during the propagation. 

In the following we use numerical simulations to 
demonstrate the effects of control light beam size on var- 
ious stationary transverse signal modes with negative de- 
tuning. The transverse spatial profile of the control field 
is chosen to be Gaussian, 



n(r) = n a e- r2 ^ 2a2 \ 



(7) 



The reported results will be exemplified with f^o = 
10 8 s~ 1 and a = 50 /xm. The parameters of the medium 
are set to be g 2 N - 
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22, 



for the atom den- 



sity of 10 i4 cm- ci 
length is 780 nm. 

Jo°° A^mnir)\ 2 dr - 



and the |0) — ► |2) transition wave- 
The eigenmodes are normalized as 
: const. 




FIG. 2: The lowest three transverse modes for m — 0: ipoi{r) 
(solid), ipo2{r) (dash), and 1^03 {r) (dot); with the probe signal 
detuning 8 — — 10 6 s _1 . See text for the other parameters. 

Figure [21 shows the radial profiles of the lowest three 
transverse modes, ijjoi(r), ip02(r), and ip03( r )i with m = 
and 8 = — 10 6 s _1 . The lowest mode tpoi(r) decreases 
monotonically, while the higher mode ipo n (r), with n > 1, 
oscillates and has (n — 1) nodes. Here, we must point 
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FIG. 3: The eigenvalues (upper panel) and the corresponding 
group velocities (lower panel), for the transverse eigenmodes 
ipoi(r) (solid lines), tpo2(r) (dashed lines), and i^03(r) (dotted 
lines), as the functions of detuning S. Other parameters are 
the same as in Fig.|5| 



V g > V g >"■ and V a >V 2n > ■■■ (8) 

at a given negative detuning, while each individual 
V g mn {8) decreases as the negative detuning reduces. 
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out that numerical calculations show that Gaussian func- 
tion can well approximate the the ground mode (4>oi(r)) 
plf . Figure depicts the eigenvalues /3 (upper panel) 
and the corresponding group velocities V g (lower panel) 
of the lowest three transverse modes as functions of de- 
tuning 5 < 0. For a given transverse mode, an increase 
in \S\ leads to an increase in (and also in wavevec- 
tor mismatch), but a decrease in the group velocity V g 
as it is calculated according to Eq. JJJJ. Note that if the 
effect of transverse spatial distribution is completely ne- 
glected, the group velocity for the present system would 
be V g = (l + g 2 N/nly 1 c = 10~ 6 c, larger than the values 
we obtained here. It is also interesting to see that at a 
given negative detuning, a higher mode is of a smaller 
group velocity. Thus, by making use of a high-order 
transverse mode of the probe signal field with a small 
negative detuning frequency, it is possible to have a slow 
light propagation in an EIT medium. 

We now study the transverse modes with nonzero or- 
bital angular momentum, i.e., the eigenmodes of Eq. (JSJ 
with m ^ 0, which represent optical vortex. Note that 
tp- m ,n = i>m,n [c-f. Eq. JSJl]. In FigJU we plot the 
"ground" and the first "excited" modes for both m = 1 
and m = 2: ipn (thin-solid), -012 (thin-dash), ip2i (thick- 
solid), and 1P22 (thick-dash), for the same EIT system 
and field parameters of Fig-El The transverse mode 
^ mt! (r) is found to have (n — 1) nodes, besides that of 
ipm=io,n(0) = 0. For a given n, i\) mn (f) with a large m 
extends to a large r. The negative detuning frequency 
dependences of these m ^ modes are presented in 
Fig-El The qualitative properties are similar as Fig. [3] 
It is found that V g xl > V 21 > V g 12 > V 22 . The group 
velocities for the transverse mode ip m n satisfy in general 



FIG. 4: Some lowest transverse modes with finite orbital an- 
gular momentum m: ipu(r) (thin-solid), ipi2(r) (thin-dash), 
ip2i{r) (thick-solid), and tp22(r) (thick-dash). Other parame- 
ters are the same as in Fig.|5| 
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FIG. 5: The eigenvalues (upper panel) and the group veloci- 
ties (lower paner) for the eigenmodes of Fig.0]as the functions 
of detuning frequency 8. 

All these calculated transverse mode velocities are 
smaller than the value of 10~~ 6 c, the group velocity with 
no consideration of the finite transverse distribution. Our 
calculations also conclude that the more focused (smaller 
a) the control field is, the smaller group velocity will be. 
On the other hand, as the size of the control field in- 
creases, the group velocity of each transverse eigenmode 
approaches to the asymptotic, transverse-effect- free value 
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of V g = (1 + g 2 N/Q.l)- 1 c, which is 10~ 6 c in the present 
EIT system of study. It also suggests that in order to 
observe the transverse effects experimentally, should not 
only the propagation length be long enough, but also the 
control held be well focused. 

It is noticed that the transverse profile of the control 
field, O(r) may vary in a realistic propagation in the 
z-direction, but it is assumed to be stationary in our 
theoretical treatment. The justification here is the fact 
that the adiabatic approximation is applicable if £l(r) is a 
slowly varying function of the propagating distance. As 
a result, the probe signal field prepared initially in an 
eigenmode ip m ,n(r) can stay in this mode at any distance 
z, as long as the adiabatic approximation is applicable. 
For a Gaussian control beam, the adiabatic approxima- 
tion requires the propagating distance smaller than the 
Rayleigh length of the control beam. In our example, 
the Gaussian beam with a spread of a — 50 /im leads to 
a Rayleigh length of 10 cm. Both the beam focusing size 
and propagation length parameters demonstrated here 
are well accessible in current experiments. Moreover, the 
stationary behavior of an eigenmode propagation of the 
signal field can sustain over a much longer distance, if 
a non-diffracting Bessel rather than a Gaussian control 
beam is used |23j ]. 

In conclusion, we have studied the transverse effects of 
the slow light propagation in EIT mediums. In terms of 



eigenequation, a criterion of localized transverse modes is 
given. Using this criterion, a complete set of transverse 
modes can be numerically calculated. Different trans- 
verse modes propagate with different group velocities; all 
are smaller than the limiting value calculated from the 
previous theory without the consideration of the trans- 
verse spatial effect. Increasing the size of the control field 
will decrease this velocity difference. Higher order, or 
larger orbital angular momentum mode will have smaller 
velocity. Our results will also play important roles in 
the coherent lig ht p ropagation control in general. For 
example, in Ref. I'2ll a transverse light guiding technique 
was used to obtain three-dimensional confinement of light 
pulses. A Gaussian approximation to the lowest trans- 
verse mode with zero orbital angular momentum (^>oi) 
was given. With the method present in this work, we 
should be able to identify the relevant transverse eigen- 
modes there. By carefully manipulating the transverse 
profile of the signal pulse, generalized and efficient real- 
izations of the controlled location and storage of photonic 
pulses will become possible. 
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